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1  Introduction 


Since  the  discovery  of  the  Toms  Effect  in  1948,  many  investigators  have  explored  the 
nature  and  causes  of  polymer-induced  friction  reduction,  using  both  experimental  and 
theoretical  methods. 

Experimental  investigations  have  focused  on  four  related  but  distinct  cases:  (1)  internal 
friction  reduction  produced  by  a  dilute  polymer  solution  of  unifonn  concentration 
flowing  through  a  pipe;  (2)  internal  friction  reduction  produced  by  a  uniform 
concentration  of  polymer  between  two  rotating  cylinders  (Couette  flow);  (3)  external 
friction  reduction  produced  by  a  unifonn  concentration  of  polymer  (‘polymer  ocean’) 
flowing  along  a  flat  plate  or  body  of  revolution;  (4)  external  friction  reduction  produced 
by  a  non-unifonn  concentration  of  polymer  ejected  into  the  flow  at  defined  locations 
along  a  flat  plate  or  body  of  revolution. 

Theoretical  investigations  fall  into  two  general  categories:  (1)  phenomenological  models 
that  describe  the  essential  features  of  the  experimental  results  in  broad  general  tenns;  (2) 
physics  models  that  attempt  to  explain  the  friction  reduction  effect  using  fundamental 
principles  of  hydrodynamics  and  the  physical  characteristics  of  the  polymers.  The 
physics  models  include  both  analytical  models  and  numerical  simulations. 

To  date,  the  theoretical  investigations  are,  at  best,  only  partly  successful.  Some  of  the 
phenomenological  models  accurately  capture  some  aspects  of  the  experimental 
observations.  The  numerical  simulations  also  succeed  in  replicating  some  experimental 
measurements.  However,  there  is  still  no  comprehensive  understanding  of  the 
mechanisms  of  polymer-induced  friction  reduction,  and,  in  more  practical  terms,  there  are 
no  good  models  that  capture  enough  detail  to  be  useful  for  engineering  design  of 
functional  drag-reduction  systems  employing  slot-ejected  polymer  where  the  polymer 
concentration  varies  substantially  downstream  from  the  slot. 

External  friction  reduction  for  use  on  submarines  and  other  marine  vehicles  is  an  area  of 
great  practical  importance.  In  such  applications,  the  polymer  must  be  ejected  at  one  or 
more  locations  along  the  vehicle’s  exterior  surface,  so  the  effects  of  non-uniform  polymer 
concentration  are  paramount.  Consequently,  a  model  that  can  describe  and  predict 
friction  reduction  produced  by  slot-ejected  polymer  along  an  external  surface  would  be 
very  useful  for  system  engineering. 

This  paper  describes  a  preliminary  version  of  such  a  model.  This  model  has  been  created 
by  assembling  partial  models  published  by  other  investigators  and  by  exploiting 
published  data  to  infer  additional  phenomenology.  Experiments  involved  cover  Professor 
Virk’s  work  published  in  1970;  the  work  of  the  Tulin  team  at  Hydronautics  through  the 
1970’s;  Vdovin  &  Smol’yakov  in  the  late  1970’s  and  early  1980’s;  Povkh  &  Pogrebnyak 
from  1979  through  1998;  Smith  and  Perkins  in  the  mid-1990’s;  Deutsch,  Petrie,  Fontaine 
&  Brungart  from  the  1990’s  to  the  present,  including  work  with  Cortana  Corporation; 
and,  most  recently,  work  by  Ceccio  and  his  team  (2006). 


1 


Results  from  this  model  are  compared  to  other  experiments  that  were  not  used  for  model 
development.  Such  comparisons  provide  limited  validation  of  the  model,  but  also 
indicate  areas  where  improvement  is  needed. 


This  model  will  be  useful  for  several  purposes. 

•  Enhanced  understanding  of  measurements  from  slot-ejected  polymer 
experiments 

•  System  design  for  slot-ejected  polymer  friction  reduction  on  marine 
vehicles 

•  Specification  of  new  experiments  to  improve  phenomenological 
component  models 

•  Identification  of  areas  for  more  focused  theoretical  investigations. 
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2  Definition  of  Concepts  and  Terminology 


Flow  is  assumed  along  a  solid  surface.  Initially  the  surface  is  treated  as  a  flat  plate  with 
finite  length,  Lp,  in  the  direction  parallel  to  the  flow.  The  surface  is  assumed  to  be  of 
infinite  extent  in  the  direction  perpendicular  to  the  flow,  making  the  problem  two- 
dimensional  in  Cartesian  geometry.  Later  versions  of  the  model  may  be  generalized  to  a 
body  of  revolution,  also  having  length,  Lp.  In  that  case,  the  problem  would  be  two- 
dimensional  in  cylindrical  geometry. 

The  flow  is  characterized  by  the  free-stream  velocity,  Lb  .  The  free-stream  velocity  is 
parallel  to  the  surface  of  the  plate  or  to  the  axis  of  revolution  of  the  body. 

The  flowing  fluid,  called  the  solvent,  is  characterized  by  its  density,  p,  and  its  kinematic 
viscosity,  v.  Alternatively,  the  dynamic  viscosity,  p  =  vp,  may  be  used  as  a  descriptor. 
Both  of  these  parameters  are  functions  of  temperature  and  salinity  of  the  solvent,  but,  to 
first  order,  are  treated  just  as  constants  of  the  material. 

In  the  vicinity  of  the  surface,  the  flow  is  characterized  by  a  boundary  layer.  At  the 
leading  edge  of  the  plate  or  body,  the  boundary  layer  may  be  a  laminar  boundary  layer. 
Some  distance  downstream  the  flow  makes  a  transition  from  laminar  to  turbulent,  and  the 
boundary  layer  becomes  a  turbulent  boundary  layer.  In  most  applications,  the  flow  of 
interest  for  purposes  of  friction  reduction  is  turbulent.  In  some  experiments  the  boundary 
layer  may  be  intentionally  tripped  by  the  experimenter  to  make  it  turbulent.  The  position, 
x,  along  the  plate  at  which  the  turbulent  boundary  layer  begins  is  called  the  virtual  origin 
of  the  boundary  layer.  The  virtual  origin  is  assigned  the  coordinate  x  =  0,  and 
downstream  positions  are  referenced  to  that  origin. 

A  common  non-dimensional  descriptor  of  the  flow  is  the  Reynolds  number.  Because  a 
length  scale  enters  into  the  Reynolds  number,  confusion  can  arise  among  Reynolds 
numbers  referenced  to  different  length  scales.  In  the  problem  to  be  considered  here,  both 
the  plate  length  and  the  longitudinal  coordinate  of  the  boundary  layer  play  roles.  The  two 
associated  Reynolds  numbers  are  distinguished  by  subscripts. 


ReL  =  UooLp/v 

(1) 

Rex  =  Uoox/v. 

(2) 

The  experimental  literature  also  includes  Reynolds  numbers  that  are  defined  on  the  basis 
of  boundary  layer  thickness  and/or  diameter  (in  the  case  of  pipe  flow).  Care  must  be 
exercised  to  distinguish  among  these  various  Reynolds  numbers. 

Within  the  turbulent  boundary  layer  there  are  distinguishable  structures.  A  viscous  sub¬ 
layer  extends  out  a  small  distance  from  the  solid  surface.  In  the  viscous  sub-layer, 
viscous  forces  are  much  larger  than  inertial  forces.  The  longitudinal  velocity  is  zero  at 
the  wall  (no-slip  condition),  and  increases,  first  linearly  across  the  viscous  sub-layer,  then 


3 


as  some  power  of  distance  from  the  wall.  In  the  outer  part  of  the  boundary  layer,  inertial 
forces  are  much  larger  than  viscous  forces,  and  the  latter  can  be  neglected.  The  velocity 
profile  in  that  region  is  a  logarithmic  function  of  distance  from  the  wall;  the  velocity 
approaches  the  free-stream  velocity  at  the  outer  edge  of  the  boundary  layer.  Between 
these  two  sub-layers  is  a  transition  or  buffer  layer  in  which  viscous  forces  and  inertial 
forces  are  comparable  in  magnitude.  In  this  report,  the  near-wall  region  is  considered  to 
include  the  viscous  sublayer  and  the  buffer  layer. 

Polymer  solutions  are  ejected  into  the  external  flow  through  one  or  more  slots  in  the  solid 
surface,  or,  less  commonly,  through  pores.  In  keeping  with  the  two-dimensional 
geometry,  the  slot  is  assumed  to  extend  indefinitely  in  the  direction  perpendicular  to  the 
flow.  Quantitative  descriptors  are  then  referenced  to  a  unit  span  of  the  slot. 

The  streamwise  coordinate  of  the  slot  is  xo,  referenced  to  the  virtual  origin  of  the 
boundary  layer.  The  downstream  distance  from  the  slot,  which  is  used  in  much  of  the 
analysis,  is  designated  as  xs.  Thus, 

X  =  x0  +  xs.  (3) 

The  performance  of  the  slot  is  characterized  by  the  ejection  angle,  <j),  and  by  the  volume 
ejection  rate  per  unit  span,  q.  An  alternate  description  is  provided  by  the  ejection 
velocity,  uej,  which  is  related  to  q  and  to  the  slot  width,  ws,  by 

q  =  wsuej.  (4) 

Since  q  has  the  dimensions  (Volume/Time)/Span  ~  Length  /Time,  it  is  often  referenced 
to  a  multiple  of  the  kinematic  viscosity  which  has  the  same  dimensions.  The  quantity 

Qs  =  67.3  v  (5) 

is  the  volume  transport  rate  per  unit  span  of  the  flow  in  the  viscous  sub-layer  of  an 
idealized  model  (from  the  wall  out  to  1 1.6  viscous  wall  units).  By  convention  it  is  often 
used  as  the  reference  for  the  volume  ejection  rate  of  the  slot;  q  is  expressed  as  a  multiple 
of  Qs. 


q  =  n  Qs.  (6) 

The  ejected  solute  (polymer)  is  characterized  by  its  initial  concentration  Co,  expressed  in 
grams/cc  or  in  weight  parts  per  million  ( wppm ).  Other  molecular  properties  of  the 
polymeric  material  are  often  used  in  physics-based  modeling  but  rarely  in 
phenomenological  modeling. 

Th q  friction  reduction,  DR,  is  the  fractional  change  in  drag  in  the  presence  of  polymer 
(Cfo  -  Cfp)  relative  to  the  drag  in  the  absence  of  polymer  (Cfo). 


DR  =  (1  -  Cfp/Cfo). 


(7) 
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Contemporary  measurement  technology  allows  measurement  of  drag  over  a  small 
localized  area.  Thus,  friction  reduction  can  be  a  local  variable 

DR  =  DR(xs),  (8) 

measured  with  respect  to  the  ejection  slot.  When  the  local  value  is  known,  an  integrated 
effective  friction  reduction  can  be  calculated. 

Finally,  much  of  the  literature  on  friction  reduction  by  slot-ejected  polymers  uses  a 
dimensionless  group  called  the  K factor.  This  factor  is  defined  as 


K  =  qc/pUooXs.  (9) 

This  dimensionless  group  can  be  factored  into  other  dimensionless  groups  in  a  variety  of 
ways.  One  factorization  emphasizes  mass  flux  and  Reynolds  number: 

K  =  (qc/p)(l/Rexs)  (10) 

Another  factorization  emphasizes  ratios  of  dimensionally  similar  fundamental  quantities. 

K  =  (c/p)(uej/Uoo)(ws/xs)  (11) 

Measured  values  of  friction  reduction  are  often  expressed  as  functions  of  K.  One  of  the 
goals  of  phenomenological  modeling  is  to  fully  explain  the  observed  relationship 
between  DR(xs)  and  K(xs). 
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3  Prior  Work 


Two  key  papers  on  slot-ejected  polymers  were  published  by  Vdovin  &  Smol’yakov  in 
1978  and  1981. 

The  central  thesis  of  these  papers  is  that  diffusion  alters  the  mean  wall  concentration  of 
polymer  as  a  function  of  distance  from  the  slot  in  repeatable  and  statistically  predictable 
ways.  They  introduced  the  concept  of  diffusion  length,  L,  defined  as  the  downstream 
distance  at  which  the  wall  concentration  drops  to  1/e  of  its  initial  value.  For  their 
experiments  free  stream  velocities  ranged  between  2.5  and  10  m/s.  Ejection  flow  rates 
ranged  between  0.77  and  15.4  cm2/s  (1.14  Qs  to  22.88  Qs). 

Vdovin  &  Smol’yakov  established  several  properties  of  the  diffusion  length. 

•  When  the  ejection  angle  is  small  (  <7°),  L  is  independent  of  the  free  stream 
velocity.  It  was  found  to  depend  only  on  the  dimensionless  group  qco/|u. 
For  small  values  of  qco/|u  (<  1 0),  L  is  proportional  to  qco/p.  For  larger 
values  of  qco/p.,  L  approaches  a  maximum  value  of  ~  70  cm. 

•  When  the  ejection  angle  is  larger  (~20°),  L  does  depend  on  the  free  stream 
velocity.  At  an  ejection  angle  of  20  experiments  on  several  different 
polymers  showed  that 


L  =  k  (qco/pUoo)  (12) 

over  the  full  range  of  measurements.  The  parameter  k  was  found  to  be 
characteristic  of  the  polymer  type,  with  kref  =  6.25  x  106  being  the 
reference  value  for  PEO  WSR-301.  For  other  polymers,  k  is  expressed  as 
k  =  k0kref.  Although  no  limiting  value  was  measured  up  to  L  ~  200  cm, 
they  did  not  rule  out  the  possibility  of  L  reaching  an  asymptotic  value  for 
very  large  values  of  qco/pUoo. 

These  observations  of  Vdovin  &  Smol’yakov  suggest  that  the  ejection  angle  plays  a  key 
role  in  the  downstream  behavior  of  the  polymer.  For  small  ejection  angles,  the  polymer 
stays  close  to  the  surface  so  that  its  behavior  is  governed  by  inner  variables  q,  Co,  p,  v. 
For  larger  ejection  angles,  the  polymer  extends  further  into  the  boundary  layer,  or  may 
actually  penetrate  the  boundary  layer,  so  that  the  outer  variable  Uoo  becomes  important  as 
well. 

As  discussed  by  Vdovin  &  Smol’yakov,  the  diffusion  length  is  to  be  interpreted  in  the 
context  of  three  ‘zones’  downstream  of  the  ejection  slot.  The  initial  zone  is  very  close  to 
the  slot,  the  solute  is  confined  to  the  near-wall  region,  and  the  wall  concentration  remains 
close  to  the  initial  concentration.  Typically  the  initial  zone  is  very  short.  In  the 
intermediate  zone,  the  diffusion  process  transports  the  solute  away  from  the  wall  until  it 
occupies  much  of  the  turbulent  boundary  layer.  The  concentration  at  the  wall  decreases 
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exponentially  in  this  zone.  Finally,  in  the  far  zone,  as  the  boundary  layer  expands 
downstream,  diffusion  continues  to  reduce  the  wall  concentration,  but  the  rate  is  algebraic 
rather  than  exponential. 

Vdovin  &  Smol’yakov  provide  analytical  models  for  the  wall  concentration  in  the 
intennediate  and  far  zones.  Those  models  are  based  on  fitting  a  large  amount  of 
experimental  data. 

In  the  intermediate  zone 


cw/co  ~  exp(-ax/Z  -  P),  (13) 

where  a  is  a  parameter  that  depends  on  the  angle  at  which  the  ejectant  is  introduced  into 
the  flow,  and  p  =  1  -  a,  which  ensures  that  cw/co=  1/e  when  x  =  L.  Only  two  ejection 
angles  were  tested  by  the  authors.  They  found  that  for  (j)  =  7°,  a  =  0.7  and  for  (j)  =  20°,  a  = 
1.8.  Note  that  cw  depends  on  Co  both  explicitly  and,  through  L,  implicitly. 

In  the  far  zone  their  data  fit  an  algebraic  expression  of  the  fonn 


cw/p  =4.16  10'5  (x/Z)'1'38.  (14) 

[Note:  Vdovin  &  SmoTyakov  (1981)  contains  an  apparent  typographical  error.  Here  10'5 
has  been  substituted  for  their  10'2  to  make  the  equation  agree  with  Figure  5  of  that  paper.] 
In  this  case  cw  does  not  depend  explicitly  on  the  initial  concentration,  but  it  depends 
implicitly  on  q,  Co,  and  Uoo  through  L. 
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4  Model  Development 

The  model  presented  herein  incorporates  three  primary  features. 

•  The  relationship  between  wall  concentration  of  the  polymer  and  the 
distance  from  the  slot.  This  relationship  is  taken  directly  from  Vdovin  and 
Smol’yakov  (1981),  with  allowance  for  blending  the  appropriate  functions 
between  the  different  diffusion  zones. 

•  The  relationship  between  wall  concentration  and  friction  reduction.  For 
small  concentration  values,  an  empirical  relationship  is  developed  in  this 
paper  on  the  basis  of  data  published  by  Petrie  et  al  (2003).  For  large 
concentration  values,  an  approximate  analytical  formula  is  based  on 
observations  of  Wu  &  Tulin  (1972),  Wu,  Fruman,  &  Tulin  (1977)  and 
Povkh  et  al  (1980). 

•  The  maximum  friction  reduction  as  a  function  of  Reynolds  number.  For 
this  calculation  the  local  friction  coefficient  for  pure  solvent  Cfo  was  taken 
from  the  empirical  function  cited  by  White  (1991).  That  function  was 
found  to  agree  very  closely  with  a  different  empirical  function  cited  by 
Schlichting  (1960).  The  local  friction  coefficient  for  maximum  friction 
reduction  Cfp  was  taken  from  Winkel  et  al  (2006)  citing  Virk  et  al  (1970). 

These  three  basic  features  were  supplemented  by  simple  subsidiary  calculations:  the 
relationship  between  K  and  xJL,  taking  account  of  possible  saturation  of  L;  the 
relationship  between  K  and  Rex,  taking  account  of  the  distance  between  the  slot  and  the 
virtual  origin  of  the  boundary  layer. 

These  components  of  the  model  are  described  in  the  next  sections. 

4.1  Wall  Concentration  vs.  x/L 

In  the  initial  zone,  we  assume 


cw  =  c0.  (15) 

In  the  intermediate  zone  we  use  the  fonnula  of  Vdovin  &  Smol’yakov 

cw  =  Co  exp(-ax/Z  -  P).  (16) 

In  the  far  zone  we  use  the  formula  of  Vdovin  &  Smol’yakov 

cw  =  p  4.16  10'5  (x/Z)'1'38.  (17) 
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To  blend  these  fonnulas  for  the  different  zones,  we  multiply  by  functions  that  roll-off 
sharply.  The  parameter  XOL1  =  .5  is  introduced  as  an  operator-selectable  parameter. 
The  roll-off  function  is  exp(-((x/Z)/XOLl)2). 

Then,  in  the  initial  zone 

ci  =  Co  exp(-((x/Z)/XOLl)2).  (18) 

In  the  intermediate  zone 

C2  =  Co  exp(-ax/Z  -  P)(l  -  exp(-((x/Z)/XOLl)2).  (19) 

In  the  far  zone 

c3  =  p  4.16  1CT5  (x/T)'L38(1  -  exp(-((x/Z)/XOL  1  )2) .  (20) 


Each  of  these  functions  has  the  proper  roll-off  behavior,  so  that  the  sum 


Cw=Ci  +  C2  +  C3  (21) 

provides  a  smooth  transition  between  zones  and  reflects  the  desired  behavior  in  each 
zone. 

4.2  Friction  Reduction  vs.  Wall  Concentration 

The  dependence  of  friction  reduction  on  the  wall  concentration  is  an  important  issue. 

The  model  depends  in  a  sensitive  way  on  the  functional  relationship  between  the  two 
variables.  Data  to  establish  such  a  relationship  is  most  productively  obtained  from  pipe 
flow  experiments  or  from  polymer  ocean  experiments  in  which  the  concentration  is 
uniform  and  presumably  known  with  reasonable  accuracy.  One  must  then  assume  that 
the  relationship  continues  to  hold  in  slot  ejection  experiments  where  the  wall 
concentration  varies  with  xs  and  other  parameters. 

Vdovin  and  SmoTyakov  present  a  functional  fonn  for  pipe  flow. 

DR  =  2.49  +  0.35  logi0(k0c/p).  (22) 

However,  they  qualify  this  expression  by  noting  that  it  holds  only  over  limited  range  of 
the  product  koc/p  (2  x  10"7  to  4  x  10'6),  essentially  equivalent  to  a  concentration  less  than 
4  wppm.  They  do  not  note  any  dependence  on  Reynolds  number  or  other  parameters. 

In  practice  it  is  well  known  that  there  is  a  maximum  friction  reduction  that  is  a  function 
of  Reynolds  number.  The  actual  friction  reduction  approaches  that  limit  rather  quickly  as 
the  polymer  concentration  increases  from  zero.  The  maximum  friction  reduction  may  be 
achieved  over  a  finite  range  of  concentration  values,  but,  as  the  concentration  continues 
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to  increase,  the  friction  reduction  again  decreases.  Any  realistic  model  should 
incorporate  those  features. 

4.2.1  Low  Concentration  Limit 

A  recent  paper  by  Petrie  et  al  (2003)  contains  some  useful  data  on  the  functional 
relationship  between  concentration  and  friction  reduction  at  the  lower  concentration 
values.  Although  the  main  thrust  of  that  paper  is  on  the  effects  of  surface  roughness,  it 
does  present  baseline  data  for  the  flow  of  a  polymer  ocean  across  a  smooth  flat  plate. 

The  data  include  variations  in  both  polymer  concentration  and  Reynolds  number.  Data 
extracted  from  the  “Figure  5”  paper  of  Petrie  et  al  are  shown  in  the  following  graph,  refer 
to  Figure  1. 


Key  features  of  Figure  1  are  that  (a)  the  friction  reduction  approaches  saturation  rapidly, 
essentially  achieving  maximum  friction  reduction  for  c  ~  20  wppm;  (b)  the  rate  at  which 
friction  reduction  approaches  its  asymptote  depends  on  the  Reynolds  number,  with  the 
largest  Reynolds  numbers  corresponding  to  the  most  rapid  growth  in  friction  reduction. 

It  would  be  tempting  to  model  the  approach  to  saturation  as  a  simple  exponential: 

DR(c)  =  DRmax(  1  -  exp(-yc)),  (23) 
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where  y  is  a  constant  that  may  depend  on  the  Reynolds  number.  Closer  examination 
shows  that  this  is  not  necessarily  a  good  model.  The  friction  reduction  rises  very  rapidly 
for  small  values  of  c,  and  then  approaches  the  asymptote  more  slowly  than  would  a 
simple  exponential.  A  more  reasonable  model  that  has  the  desired  properties  is 

DR(c)  =  DRmax(  1  -  exp(-y(Re)cpp)),  (24) 

where  pp  is  a  constant  of  less  than  one  and  y  appears  to  depend  on  flow  velocity. 

Least  squares  fits  to  the  data  of  Petrie  et  al  for  each  Reynolds  number  shows  that  0.4  < 
PP  <  0.6,  with  no  systematic  dependence  of  PP  on  Reynolds  number.  The  mean  value 
over  the  range  of  Reynolds  numbers  is  <PP>  =  0.54.  The  coefficient  y  does  depend  on 
the  flow  velocity  in  a  systematic  way.  A  least  squares  fit  for  that  parameter  gives 

y  =  1.05  (Reref/Re)'83.  (25) 

Consequently,  the  scaling  c°'54/Re0'83  should  collapse  the  data.  Figure  2  shows  this 
scaling.  The  function  1  -  exp(-1.05(Reref/Re)°'83c0'54)  is  shown  for  comparison. 


Figure  2:  Petrie,  et  al.  scaled  data  for  various  Reynolds  numbers:  4e6  <  Re  <  l.le7 
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An  important  aspect  of  the  data  of  Petrie  et  al  is  that  they  used  ReL,  not  Rex  for  the 
Reynolds  number.  Thus,  the  only  variable  factor  in  their  Reynolds  number  is  the 
velocity,  so  the  friction  reduction  scaling  could  be  written  instead  as 

(Urei/Uoo)0-83c0-54.  (26) 

This  fonn  emphasizes  that  concentration  and  velocity  scale  with  different  exponents.  As 
a  consequence,  the  K  factor,  in  which  c  and  1/U®  appear  with  the  same  exponent,  cannot 
be  expected  to  completely  collapse  the  data. 

Extraction  of  scaling  exponents  from  the  data  of  Petrie  et  al  is  rather  crude,  and  the  exact 
values  of  the  exponents  may  vary.  Clearly,  a  much  larger  set  of  polymer  ocean  data 
would  be  needed  to  arrive  at  a  fully  reliable  estimate. 

Two  sets  of  pipe  flow  data  support  an  estimate  of  PP  ~  .5.  Unfortunately,  those  data  sets 
do  not  contain  information  on  the  concurrent  effects  of  varying  concentration  and 
velocity  (or  Reynolds  number).  Virk  &  Baher  (1970)  found  an  effect  proportional  to  c'/2 
over  two  orders  of  magnitude  in  the  variation  of  c.  Data  from  Warholic  et  al  (1999)  is 
shown  in  Figure  3.  The  best  fit  to  that  data  at  Re  ~  20,000  gives  pp  =  0.527. 
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In  the  present  model,  the  friction  reduction  will  be  incorporated  as 

DR  =  DRmax  ( 1  -  exp(-y(Uref/U*)scpp)),  (27) 

with  nominal  values  y  =  1.05,  8  =  0.83,  pp  =  0.54,  and  Uref  =  457  cm/sec  estimated  from 
Petrie  et  al.  (pp  is  used  here  to  distinguish  this  parameter  from  the  p  used  by  Vdovin  & 
Smol’yakov.)  DRmaxis  a  function  of  Reynolds  number,  as  discussed  in  Section  4.3. 

4.2.2  High  Concentration  Limit 

Wu  &  Tulin  (1972)  first  showed  that  friction  reduction  performance  decreases  for  high 
concentrations  of  polymer.  For  a  polymer  ocean  flowing  along  a  flat  plate,  they  showed 
that  maximum  friction  reduction  occurred  for  c  ~  1 00  wppm.  For  slot  ejection  into  a 
flow  with  Uoo  =  8  ft/sec,  they  showed  that  integrated  friction  reduction  for  the  whole  plate 
increased  with  concentration  up  to  c  ~  500  wppm,  then  remained  constant  up  to  the  limits 
of  their  experiment  at  c  =  1000  wppm.  Subsequent  slot-ejection  work  (Wu,  Fruman  & 
Tulin,  1977)  extended  the  concentration  range  to  5000  wppm.  That  work  showed 
decreased  friction  reduction  at  the  end  of  the  approximately  3-meter  plate  for  1000,  3000, 
and  5000  wppm,  compared  to  best  performance  in  the  range  250  -  500  wppm. 

Povkh  et  al  (1979)  used  small  diameter  tubes  to  measure  friction  reduction  for  various 
concentrations  and  various  molecular  weights  of  PEO  WSR  301.  For  the  molecular 
weight  M  =  106,  roughly  comparable  to  that  of  the  polymer  used  by  Wu  &  Tulin,  they 
also  found  the  maximum  friction  reduction  to  occur  at  c  ~  100  wppm.  However,  the 
decrease  in  friction  reduction  at  higher  concentrations  was  somewhat  slower  than  that 
reported  by  Wu  &  Tulin.  They  also  noted  a  weak  dependence  of  the  optimum 
concentration  on  Reynolds  number. 

Both  sets  of  data  are  rather  sparse.  An  approximate  model  for  the  high  concentration 
behavior,  without  consideration  of  the  effects  of  molecular  weight  or  Reynolds  number 
is: 


DR  =  DRmax  (croll/(c  +  Croll)),  (28) 

where  the  roll-off  concentration,  croii,  is  an  empirical  parameter.  A  value  of  cron  ~  3500 
wppm  is  used  in  the  present  illustration  of  the  model.  That  value  gives  a  reasonable  fit  to 
the  data  of  Povkh  et  al,  but  overestimates  the  data  of  Wu  &  Tulin.  In  the  implementation 
of  the  model,  cron  is  an  operator-selectable  parameter. 

Combination  of  the  low-concentration  and  high-concentration  models  gives  the  final 
version  used  in  the  present  calculations: 

DR  =  DRmax  ( 1  -  exp(-y(UrekU00)5cpp))(croii/(c  +  cron)).  (29) 
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4.3  Maximum  Friction  reduction 


Many  experiments  in  different  experimental  configurations  have  shown  that  there  is  a 
limiting  (minimum)  friction  factor,  Cfp,  that  can  be  achieved  by  the  use  of  polymer 
solutions.  That  limit  is  a  function  of  Reynolds  number.  Similarly,  the  friction  factor  in 
the  absence  of  polymer,  Cfo,  is  also  a  function  of  Reynolds  number. 

The  maximum  friction  reduction  is 


DRmax(Rex)  =  1  -  Cfp(Rex)/Cf0(Rex).  (30) 

The  friction  factors  are  available  as  empirical  functions  estimated  from  extensive 
experimental  data.  Cfo  is  given  as  an  explicit  function  by  White  (1991) 

Cfo  =  0.455/ (loge(0 . 06Rex))2 .  (31) 

An  alternative  empirical  function  given  by  Schlichting  (1960)  produces  a  nearly  identical 
numerical  result  over  the  range  of  Reynolds  numbers  of  interest. 

Cfo  =  (2  log10(Rex)-0.65)'2'3.  (32) 

The  expression  from  White  is  used  in  the  present  model. 

Cfp  is  given  implicitly  by  a  function  cited  by  Winkel  et  al  (2006)  referring  to  the  work  of 
Virk  et  al  (1970). 

Qp-1/2  =  19  log10(Cfp  Rex)  -  38.1  (33) 

In  this  model  Cfp  is  calculated  as  a  function  of  Rex  by  fixed  point  iteration  of  this 
equation.  Note  that  this  is  different  from  the  conventional  formula  for  the  Virk 
asymptote  in  pipe  flow  which  is  based  on  the  diametral  Reynolds  number  Rcd,  not  Rex. 

It  is  also  important  to  note  that  the  empirical  Virk  asymptote  is  based  on  studies  of  shear 
flow.  Topologically,  turbulent  shear  flow  is  -50%  extensional  flow  and  ~50%  rotational 
flow.  A  different  type  of  flow,  in  which  the  proportions  of  extensional  and  rotational 
flows  differ,  might  have  a  different  maximum  friction  reduction  asymptote. 

4.4  Other  Functions 

Two  final  relationships  are  needed  to  close  the  model.  It  is  desired  to  express  the  results 
in  terms  of  the  K  factor.  The  wall  concentration  is  expressed  in  tenns  of  the  ratio  xs/L  in 
the  formulas  of  Vdovin  &  Smol’yakov.  The  maximum  friction  reduction  is  expressed  in 
tenns  of  the  Reynolds  number  Rex,  based  on  distance  from  the  virtual  origin  of  the 
boundary  layer.  It  is  necessary  to  relate  K  to  Rex  and  xJL. 
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4.4.1  K  vs  xs/L 


The  relationship  between  K  and  xs/X  recognizes  the  possibility  that  X  may  reach  some 
upper  limit  Xmax.  By  definition 


K  =  qc0/pUooXs  (34) 

Multiplication  by  1  =  (X/X)(Xmax/Xmax)  gives  K  =  (X/xs)  ( Xmax/X )  (qc0/pU,Xmax). 

When  qco/pUoo  is  small,  X  is  proportional:  X  =  k  qco/pUoo.  If  X  approaches  an  asymptotic 
value  Xmax  when  qco/pUoo  becomes  large,  X  can  be  modeled  as 

L  —  f  max  (1  —  exp(-k  qCo/pUccZ'max))-  (35) 


Then  xs/X  can  be  written 


xs/X  =  (1/kK)  F,  (36) 

where  the  correction  factor  F  is 

F  =  (kqco/pUoofmax)/  ( 1  -  cxp(-kqco/pU,Xmax)).  (37) 

When  kqco/pUooXmaxis  small,  F  =  1,  xs/L  =  (1/kK),  and  no  correction  is  necessary.  When 
k q c 0/ pU  /Xmax  becomes  large,  F  ->  kqco/pUooXmax  and  additional  factors  of  q,  Co,  and  1/Uoo 
are  introduced  into  xJL. 

4.4.2  K  vs.  Rex 

The  principal  distinction  here  is  the  choice  of  origins  for  K  and  Rex.  K  is  referenced  to 
the  position  of  the  slot  and  Rex  is  referenced  to  the  virtual  origin  of  the  turbulent 
boundary  layer.  Since  x  =  xo  +  xs,  multiplication  by  U  ,/v  gives 

Rex  =  Rexo  +  (qc0/|u)(l/K),  (38) 

where  Rexo  =  UooXo/v  is  the  Reynolds  number  based  on  the  distance  from  the  virtual  origin 
of  the  boundary  layer  to  the  slot.  This  number  will  obviously  change  with  flow  velocity, 
and,  depending  on  the  location  of  the  slot,  can  become  quite  significant. 

4.5  Documentation 

The  complete  model  is  coded  in  MatLab.  In  the  model,  all  units  are  cgs,  so  distances  are 
in  centimeters  and  velocities  are  in  centimeters  per  seconds. 
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5  Behavior  of  the  Model 


The  model  requires  a  large  number  of  inputs.  These  can  be  grouped  into  three  categories: 
(1)  established  physical  parameters:  p,  v,  p  that  characterize  the  solvent;  (2)  empirical 
parameters  based  on  prior  experiments:  k0,  kref,  Zmax,  Uref,  the  reference  velocity  in 
Petrie’s  experiments,  the  parameters  used  to  fit  the  data  of  Petrie  et  al,  y,  8,  (3(3,  and  cron, 
the  parameter  used  to  fit  the  high  concentration  behavior;  (3)  experimental  parameters 
that  characterize  a  specific  experiment:  free  stream  velocity  Ux,  initial  concentration  Co, 
ejection  angle  parameter  a,  ejection  rate  q,  and  position  of  the  slot  with  respect  to  the 
virtual  origin  of  the  turbulent  boundary  layer  xo.  The  first  group  can  be  taken  as  fixed 
parameters  that  do  not  affect  the  functioning  of  the  model.  The  second  group  needs  to  be 
investigated  to  determine  the  sensitivity  of  the  model  to  the  empirical  parameters. 
Although  the  empirical  parameters  have  been  fit  to  prior  experimental  data,  it  is  still 
necessary  to  evaluate  them  and  ensure  that  the  results  of  the  model  are  consistent  with 
new  experiments.  This  process  can  be  viewed  as  a  form  of  validation.  Finally,  the  model 
can  be  used  to  investigate  variations  in  the  experimental  conditions  expressed  by  the  third 
group  of  parameters. 


5.1  Sensitivity  to  Empirical  Parameters 

In  this  section,  we  choose  a  nominal  set  of  parameters  and  operating  conditions  as  a 
baseline,  and  then  investigate  the  sensitivity  of  the  model  to  the  empirical  parameters, 
one-by-one.  The  joint  sensitivities  of  pairs  of  parameters  will  be  considered  in  the  future. 

The  nominal  operating  conditions  are: 

Ua)  =  900  cm/sec; 

Co  =  1000  wppm; 
q  =  5  Qs; 

a  =  1.8  (slot  angle  =  20°); 
x0  =  70  cm. 

The  nominal  empirical  parameters  are 

k0=  1  (PEO  WSR-301); 
kref  =  6.25  x  106; 

Lmax  270  cm, 

Uref  =  457  cm/sec; 
y  =  0.80; 

8  =  0.83; 

(3(3  =  0.54; 
cron  =  3500  wppm 
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Figure  4  shows  the  friction  reduction  as  a  function  of  K  for  the  nominal  case. 


Figure  4:  Nominal  baseline  for  sensitivity  analysis 


Figure  4  is  typical  of  experimental  data  from  slot-ejected  polymer  experiments.  The 
friction  reduction  is  insignificant  for  K  <  10"9.  It  increases  steadily  and  is  more  or  less 
proportional  to  logio(K)  in  the  range  10"  <  K  <  10"  .  It  then  reaches  a  plateau  and  may 
decline  in  the  near-slot  region. 
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The  factors  that  contribute  to  this  shape  are  shown  in  Figure  5.  Contributions  from  the 
initial  zone,  the  intermediate  zone,  and  the  far  zone  are  compared  with  the  maximum 
friction  reduction  predicted  by  Virk’s  asymptote.  For  this  particular  set  of  operating 
conditions,  the  contribution  from  the  far  zone  is  equal  to  the  contribution  from  the 
intermediate  zone  at  K  ~  6  x  10'  .  For  smaller  values  of  K,  the  far  zone  dominates.  For 
larger  values  of  K,  where  the  contributions  from  the  intermediate  and  initial  zones  might 
be  expected  to  dominate,  Virk’s  asymptote  is  a  governing  factor.  In  this  specific 
example,  the  intermediate  zone  has  a  significant  influence  only  over  a  very  limited  range 
of  K.  In  the  initial  zone,  the  influence  of  the  decrease  in  performance  due  to  high 
polymer  concentration  is  apparent. 


Figure  5:  Components  of  total  friction  reduction 
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Figure  6  shows  the  total  concentration  at  the  wall  as  a  function  of  K  and  compares  the 
total  with  the  model  results  for  the  three  zones.  For  this  example,  the  cross-over  point 

o 

between  the  intermediate  zone  and  the  far  zone  at  K  ~  6.5  10'  is  obvious.  The  cross-over 
point  between  the  initial  zone  and  the  intermediate  zone  occurs  around  K  ~  4.5  10'7. 


Figure  6:  Components  of  wall  concentration 


The  following  sections  examine  the  sensitivity  of  the  model  to  each  of  the  parameters. 
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5.1.1  Sensitivity  to  k0 


The  ko  parameter  was  introduced  by  Vdovin  &  SmoTyakov  to  compare  the  performance 
of  different  polymer  types.  By  definition  ko  =  1  refers  to  PEO  WSR-301.  Any  polymer 
with  poorer  performance  will  have  a  value  of  ko  <  1 .  Any  polymer  with  superior 
performance  will  have  a  value  of  ko  >  1.  Note  that  Vdovin  &  SmoTyakov  relate  this 
standard  value  to  ‘fresh’  polymer,  and  they  provide  data  that  indicate  that  ‘aged’  WSR- 
301  can  have  a  ko  value  of  ~0.2.  As  Figure  7  shows,  performance  of  the  model  is 
sensitive  to  ko.  The  general  effect  is  to  shift  the  curves  left  or  right  without  significant 
change  of  shape. 
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5.1.2  Sensitivity  to  Lmax 


The  parameter  Lmax  was  introduced  to  account  for  the  possibility  that  the  diffusion  length, 
L,  saturates  or  reaches  an  asymptotic  value  as  other  experimental  parameters  change.  As 
Vdovin  &  Smol’yakov  showed,  Lmax  ~  70  cm  when  the  ejection  angle  is  small.  However, 
when  the  ejection  angle  was  large,  they  did  not  find  any  limiting  behavior  of  L.  Figure  8 
shows  that  the  model  behavior  is  not  sensitive  to  Lmax  unless  Lmax  assumes  unrealistically 
small  values.  In  practical  terms,  this  implies  that 

xJL  ~  1/kK  (39) 

almost  everywhere  without  the  need  for  the  correction  factor  F  of  Equation  (37). 
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5.1.3  Sensitivity  to  Uref 


Uref  is  a  parameter  taken  directly  from  Petrie  et  al,  and  there  is  no  reason  to  assume  that  it 
is  in  error.  As  Figure  9  shows,  the  model  is  not  sensitive  to  25%  -  35%  errors  in  Uref. 

The  only  effect  is  a  slight  shift  of  the  curves  left  or  right  without  change  of  shape.  Thus 
the  nominal  value  of  this  parameter  appears  to  be  adequate  for  the  model. 


22 


5.1.4  Sensitivity  to  y 


The  parameter  y,  which  appears  in  the  exponent  that  governs  the  approach  to  saturation  as 
a  function  of  wall  concentration,  has  a  modest  influence  of  the  model.  In  general, 
increasing  y  increases  the  apparent  performance  of  the  polymer.  As  Figure  10  shows, 
changes  in  y  cause  a  shift  of  the  whole  curve  left  or  right  without  significant  change  of 
shape. 
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5.1.5  Sensitivity  to  8 


The  parameter  8  is  the  power  law  exponent  that  describes  the  contribution  of  the  velocity 
in  the  approach  to  saturation  as  a  function  of  wall  concentration.  Figure  1 1  shows  that 
the  model  is  not  sensitive  to  changes  of  8  within  a  reasonable  range  around  its  nominal 
value. 
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5.1.6  Sensitivity  to  (3(3 


The  parameter  pp  is  the  power  law  exponent  of  the  concentration  in  the  approach  to 
saturation  from  below.  Figure  12  shows  that  this  parameter  has  a  major  impact  on  the 
magnitude  and  shape  of  the  DR  vs.  K  curve.  In  particular,  it  changes  the  slope  of  the 
curve  in  the  diffusion  far  zone.  This  parameter  provides  a  sensitive  test  of  the  correctness 
of  the  model.  It  is  thus  very  important  to  have  a  reliable  estimate  of  the  value  of  pp  from 
independent  sources.  The  data  of  Petrie  et  al  from  polymer  ocean  experiments  is  limited 
in  scope  and  accuracy.  Examination  of  additional  data  from  other  polymer  ocean 
experiments  and  from  pipe  experiments  should  improve  the  estimate  of  pp. 

At  a  deeper  level,  it  will  also  be  important  to  understand  why  this  scaling  parameter  is 
approximately  1/2,  and  to  justify  the  assumption  that  the  scaling  approach  to  saturation 
observed  in  polymer  ocean  and  pipe  experiments  also  applies  to  external  flow  conditions. 


Figure  12:  Sensitivity  to  Double  Beta  (PP) 
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5.1.7  Sensitivity  to  cron 


The  value  of  the  roll-off  parameter  for  high-concentration  effects,  cron,  has  a  significant 
impact  on  the  model  in  the  near-slot  region  where  the  initial  polymer  concentration  is 
high,  refer  to  Figure  13. 


5.2  Influence  of  Experimental  Conditions 

Polymer  friction  reduction  experiments  have  been  conducted  on  the  external  surface  of 
flat  plates,  bodies  of  revolution,  and  operational  marine  vehicles.  For  purposes  of 
comparing  experimental  conditions,  the  flat  plat  provides  the  simplest  example.  The 
principal  parameters  that  can  be  varied  in  the  course  of  an  experimental  investigation  are: 
U„,  the  free  stream  velocity;  xo,  the  longitudinal  coordinate  of  the  slot  with  respect  to  the 
virtual  origin  of  the  turbulent  boundary  layer;  the  slot  ejection  angle,  4>,  manifest  through 
the  angle  parameter,  a;  the  initial  concentration  of  ejected  polymer,  Co;  and  the  polymer 
ejection  rate,  q. 

The  influences  of  these  experimental  parameters  on  the  model  are  discussed  next. 

As  in  Section  5.1,  the  set  of  nominal  empirical  parameters  and  the  set  of  nominal 
operating  conditions  provide  the  baseline  case.  The  effects  of  the  experimental 
parameters  are  demonstrated  by  varying  them,  one  by  one,  from  their  nominal  values. 
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5.2.1  Influence  of  Free  Stream  Velocity 


The  influence  of  free  stream  velocity  is  shown  in  Figure  14.  Two  separate  effects  are 
apparent. 

In  the  diffusion  far  zone  region,  there  is  a  significant  separation  between  performance 
curves  at  different  values  of  Uoo.  Superior  performance  is  predicted  by  the  model  at 
smaller  values  of  Uoo.  That  separation  is  a  consequence  of  the  scaling  of  the  approach  to 
saturation,  (Ure/Uoo)5  The  fact  that  both  Uoo  and  c  obey  power  laws  with  different 
exponents  ensures  that  the  ratio  c/Ua,  does  not  collapse  to  a  common  curve. 

In  the  region  controlled  by  Virk’s  asymptote,  the  order  is  reversed,  and  superior 
performance  is  attained  at  higher  values  of  Uoo.  That  is  because  the  maximum  friction 
reduction  predicted  by  Virk’s  asymptote  is  an  increasing  function  of  Reynolds  number, 
hence  of  U®. 


Figure  14:  Variation  with  free  stream  velocity  (m/s) 


5.2.2  Influence  of  Slot  Position  re  Virtual  Origin 

The  parameter  xo  only  influences  friction  reduction  performance  in  the  region  controlled 
by  Virk’s  asymptote.  This  is  also  a  consequence  of  the  dependence  of  the  asymptote  on 
Reynolds  number.  For  a  given  value  of  K,  the  operative  Reynolds  number  in  the 
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calculation  of  the  maximum  friction  reduction,  Rex,  is  increased  by  the  Reynolds  number 
Rexo,  which  is  a  linear  function  of  x0.  Thus,  as  x0  increases  (greater  separation  between 
the  virtual  origin  and  the  slot  position),  the  maximum  friction  reduction  also  increases,  as 
illustrated  in  Figure  15. 


Figure  15:  Variation  with  Slot  Position  x0  (cm) 


5.2.3  Influence  of  Slot  Angle  Parameter  a 

The  model  does  not  yet  contain  a  relationship  between  the  ejection  angle  and  the 
parameter  a.  Vdovin  &  Smol’yakov  only  report  results  for  two  ejection  angles.  Their 
results  show  that  a  is  an  increasing  function  of  <j).  Values  of  a  for  ejection  angles  greater 
than  their  maximum  angle  of  20  must  be  extrapolated. 

Sensitivity  studies  with  the  present  model  show  that  the  predicted  perfonnance  of  the 
polymer  depends  very  strongly  on  a.  For  that  reason,  accurate  knowledge  of  a  over  a 
realistic  range  of  ejection  angles,  as  well  as  a  fundamental  understanding  of  the  physics 
of  varying  ejection  angles,  becomes  very  important. 

Figure  16  shows  that  small  values  of  a  extend  the  intermediate  diffusion  zone  to 
significantly  smaller  values  of  K.  Part  of  that  extended  range  is  controlled  by  Virk’s 
asymptote,  allowing  high  friction  reduction  and  superior  performance  over  a  greater 
distance  downstream  of  the  slot.  Where  the  performance  is  controlled  by  the  exponential 
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damping  function  of  the  intermediate  zone,  the  slope  of  the  performance  curve  is  very 
steep.  For  smaller  values  of  K,  the  high-slope  intermediate  zone  blends  into  the  gentler 
algebraic  slope  of  the  far  zone. 

The  model  implies  that  smaller  ejection  angles  and  smaller  a  values  would  show  superior 
performance  over  an  even  wider  region  of  K  space. 


Figure  16:  Variation  with  Slot  Angle  Parameter  Alpha  (a) 


5.2.4  Influence  of  Initial  Concentration  c0 

According  to  the  results  of  Vdovin  and  Smol’yakov,  the  wall  concentration  is  explicitly 
proportional  to  the  initial  concentration  in  the  initial  and  intermediate  zones,  but  not  in 
the  far  zone.  That  dependence  of  their  model  is  reflected  in  Figure  17,  which  shows 
significant  dependence  of  performance  on  Co  in  the  intermediate  zone.  The  dependence 
essentially  vanishes  in  the  far  zone,  where  any  dependence  on  Co  is  implicitly  contained 
in  L.  But  the  strongest  dependence  on  Co  occurs  in  the  initial  zone,  where  the  effects  of 
entanglement  become  manifest  for  large  values  of  concentration.  The  dependence  on  Co 
is  also  impacted  by  the  choice  of  cron.  In  the  example  shown,  cron=  3500  wppm.  This 
causes  a  large  decrease  in  friction  reduction  when  Co  =  4000  >  cron,  but  produces  no 
obvious  effect  when  Co  =  100  «  cron. 
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Figure  17:  Variation  with  Initial  Polymer  Concentration  c0  (wppm) 


5.2.5  Influence  of  Polymer  Ejection  Rate  q 

In  the  same  spirit,  the  friction  reduction  performance  depends  weakly  on  the  ejection  rate, 
q,  only  in  the  initial  and  intermediate  zones.  Figure  18  shows  that  ejection  rate  does  not 
influence  performance  in  the  far  zone.  In  that  zone  the  K  factor  collapses  the  q 
dependence  to  a  single  line  that  incorporates  the  product  qco/px,  but,  as  shown  in  Section 
5.2.1,  does  not  properly  collapse  with  l/lf*,.  As  discussed  in  Section  7  below,  q  is 
determined  by  the  ejection  velocity  and  the  slot  width. 
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Figure  18:  Variation  with  Polymer  Ejection  Rate  (q) 
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6  Comparison  with  Published  Experiments 


In  this  section,  the  model  is  exercised  to  compare  its  predictions  with  data  published  by 
various  experimental  groups.  These  data  are  independent  of  the  data  used  to  build  the 
model.  Consequently,  these  comparisons  should  provide  a  good  test  of  the  model. 

6.1  Data  of  Winkel  et  al  (2006) 

Recent  experiments  were  reported  by  Winkel  et  al  (2006).  Their  results  are  of  particular 
interest  because  they  show  a  strong  effect  of  free  stream  velocity.  Data  taken  at  different 
values  of  Uoo  do  not  collapse  onto  a  single  curve  in  the  DR  vs.  K  plot. 

The  experiments  were  conducted  in  the  William  B.  Morgan  Large  Cavitation  Channel  in 
Memphis,  TN.  Experimental  conditions  were  as  follows. 

Polymer  type: 

PEO  WSR-301,  WSR-N-60K,  WSR-308 
Slot  Position  re  Leading  Edge: 

xo:  132  centimeters 
Drag  Measurements  re  Slot  Position: 

xs:  -25,  64,  209,  611,  791,  1018  centimeters 
Free  Stream  Velocity: 

Uoo:  600,  1200,  1800  centimeters/second 
Polymer  Concentration: 

Co:  1000,  2000,  4000  wppm 
Ejection  Rate: 

q:  2Qs,  4Qs,  lOQs 
Ejection  Angle: 

(|>  =  25°. 

Extensive  amounts  of  data  were  collected,  and  several  of  the  data  sets  were  presented  in 
the  DR  vs.  K  format.  Of  particular  interest  are  (1)  comparisons  among  different  free 
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stream  velocities  for  PEO  WSR-301  (Winkel  et  a/’s  Figure  16,  reproduced  here  as  Figure 
19),  and,  (2)  comparisons  among  different  polymer  types  for  various  free  stream 
velocities. 


Figure  19:  |from  Winkel,  et  al.,  Proc.  26th  Symposium  on  Naval  Hydrodynamics,  (2006)] 


To  make  comparisons  with  the  published  data,  the  model  was  exercised  with  the  nominal 
empirical  parameters.  Those  parameters  were  then  adjusted  to  provide  the  best  fit  to  the 
data  of  Winkel  et  al.  Large  adjustments  were  found  to  be  unnecessary.  All  of  the 
experimental  parameters  were  specified  except  a.  The  largest  ejection  angle  tested  by 
Vdovin  &  Smol’yakov  was  20°,  corresponding  to  a  =  1.8.  In  the  present  case,  the  a 
value  was  extrapolated  to  -25°  by  scaling  according  to  sine(c))),  giving  a  2.25 .  At 
present  the  sine(4>)  scaling  is  a  heuristic  without  a  rigorous  basis  in  theory. 

Figure  20  compares  the  output  of  the  model  with  the  data  of  Winkel  et  al.  for  three 
different  free  stream  velocities.  The  match  appears  to  be  good.  DR  values  at  K  values  K 
=  10'  and  K  =  10'  are  accurate  for  each  free  stream  velocity.  The  irregular  spacing 
between  the  different  velocity  curves  is  also  accurately  predicted.  The  decreased 
perfonnance  due  to  high  concentrations  in  the  near-slot  region  (initial  zone)  is  also 
predicted. 
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Figure  20:  Comparison  of  Model  with  Data  of  Winkel,  et  al. 


Differences  in  polymer  type  are  expressed  through  the  empirical  parameter  ko  of  Vdovin 
&  Smol’yakov.  The  reference  value  is  ko  =  1  for  WSR-301.  Based  on  pipe  flow  data, 
WSR-308  has  a  value  ko  ~  2.  The  value  for  WSR  N60-K  is  not  known.  A  value  ko  =  0.2 
was  used  to  give  a  reasonable  match  to  the  data  of  Winkel  et  al. 

Figure  22  shows  the  model  output  for  the  three  polymer  types  at  Uoo  =  6  m/s.  The  curves 
for  301  and  308  accurately  reflect  the  data  in  Winkel  et  al  [Winkel  et  a/’s  Figure  15(A) 
reproduced  here  as  Figure  21]  without  any  parameter  adjustment.  The  curve  for  N60-K 
accurately  matches  the  data  after  assignment  of  ko  =  0.2  to  that  polymer. 
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Figure  21:  |from  Winkel,  et  al.,  Proc.  26th  Symposium  on  Naval  Hydrodynamics,  (2006)) 


Figure  22:  Variation  with  Polymer  Type:  6  m/s 
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6.2  Data  of  Cortana  Corporation 


Experiments  were  conducted  at  the  Applied  Research  Laboratories  of  Pennsylvania  State 
University,  using  an  ejection  slot  designed  by  Cortana  Corporation. 

Experimental  conditions  were  as  follows. 

Polymer  Type: 

PEO  WSR301;  WSR309 

Slot  Position  re  Virtual  Origin  (Tripped  BL): 

xo:  27.9  centimeters 

Drag  Measurements  re  Slot  Position: 

xs:  7.7,  200,  329,  474,  646  centimeters 

Free  Stream  Velocity: 

Ux,:  457,  762,  1067  centimeters/second 

Polymer  Concentration: 

Co:  500,  1000,  2000  wppm 

Ejection  Rate: 

q:  5Qs,  lOQs,  20Qs 

Ejection  Angle: 

cj)  =  <  5°. 

Comparisons  between  measurements  and  model  predictions  are  shown  in  Figure  23, 
where  the  model  is  modified  to  adhere  to  the  observation  of  Vdovin  &  Smol’yakov  that: 
for  small  ejection  angles,  the  behavior  of  the  polymer  is  governed  by  inner  variables 
alone,  not  by  outer  variables  such  as  Uyj.  That  is,  in  the  intennediate  and  far  zones,  the 
data  collapse  to  a  single  function  of  K,  independent  of  free  stream  velocity. 

An  important  feature  of  this  case  is  the  effect  of  the  small  slot  angle.  The  use  of  a  small 
value  for  the  a  parameter  extends  the  impact  of  the  intermediate  zone  toward  smaller  K 
values  and  predicts  a  steeper  slope  in  the  region  10'  <K<4xlO'.  Those  predictions 
match  the  ARL/Cortana  experimental  data  very  well.  Comparison  with  the  data  of 
Winkel  et  al  (Figure  20)  shows  that  their  large  slot  angle  (25°)  minimizes  the  role  of  the 
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intermediate  zone;  their  shallower  slope  in  the  region  10'  <  K  <  10'  indicates  behavior 
characteristic  of  the  far  zone  over  most  of  their  measurement  range. 

The  Cortana  data  exhibits  large  variability  in  the  initial  zone,  that  is,  in  the  near-slot 
region  (K  >  2  x  10'7),  including  some  performance  that  exceeds  the  predicted  Virk 
asymptote.  A  separate  analysis  in  the  following  section  indicates  that  the  ejection 
velocity  (related  to  q  and  slot  width)  has  a  significant  impact  in  the  near-slot  region. 


K 

Figure  23:  WSR  301  Cortana  Fluidics  Ejector;  Model:  Inner  Variable  Flypothesis 
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7  Near-Slot  Effects 


Discrepancies  between  the  model  and  the  Cortana  data  tend  to  be  largest  for  large  values 
of  K.  Those  K  values  correspond  to  small  values  of  x,  very  close  to  the  injection  slot.  In 
that  region,  there  is  a  large  amount  of  scatter  in  the  friction  reduction  performance; 
however,  the  scatter  is  not  totally  random,  and  appears  to  be  correlated  with  q  and  Co, 
variables  that  are  related  to  the  properties  of  the  polymer  at  the  point  of  ejection. 

In  the  Cortana  WSR  301  data  taken  on  the  GTV  by  ARL/PSU,  there  were  14  data  sets  for 
measurements  at  a  location  7.7  centimeters  from  the  slot.  Friction  reduction  varied  from 
.38  to  .83,  from  well  below  to  slightly  above  Virk’s  asymptote  for  that  location  and 
Reynolds  number.  The  free  stream  velocity  did  not  appear  to  have  any  significant 
influence  on  the  friction  reduction  for  these  data  sets.  The  ejection  rate  q  had  a  strong 
positive  influence,  and  the  initial  concentration  of  polymer,  Co,  had  a  weak  and 
statistically  insignificant  influence  on  friction  reduction. 

A  linear  regression  model  of  DR  vs  q/Qs  and  Co  in  the  near-slot  region  produced  a  fit  that 
accounted  for  66%  of  the  variance: 

DR  =  .2826  +  .0194  q/Qs  +  .000039  c0.  (40) 

q/Qs  alone  accounted  for  63%  of  the  variance.  These  calculations  show  that  q/Qs  plays  a 
dominant  role  in  the  friction  reduction  perfonnance  of  polymer  solutions  near  the  ejection 
slot,  and  raise  an  important  question  about  the  physical  mechanism  behind  this  behavior. 

At  present  an  accurate  quantitative  model  of  the  near-slot  environment  is  not  yet 
available.  Nonetheless,  there  are  enough  observations  in  the  literature  to  support 
qualitative  arguments  and  to  build  an  approximate  model. 

First,  it  is  widely  accepted  that  enhanced  drag-reduction  performance  of  a  polymer 
solution  is  observed  when  the  individual  polymer  molecules  are  in  the  extended 
configuration,  rather  than  in  the  coiled  configuration.  Thus,  the  coil-stretch  transition  is 
an  important  marker  of  friction  reduction  perfonnance. 

Quantitative  description  of  this  transition  typically  uses  a  two-state  model:  transition 
from  the  coiled  configuration  to  the  stretched  configuration  is  caused  by  elongational 
forces  in  the  fluid  flow;  reversion  to  the  coiled  state  is  caused  by  elastic  forces  of  the 
molecule  itself  when  the  extensional  fluid  forces  are  temporarily  reduced.  The  quasi¬ 
equilibrium  condition  of  a  polymer  molecule  in  a  particular  flow  is  detennined  by  the 
relative  time  scales  on  which  these  competing  processes  operate.  An  extended  polymer 
molecule  relaxes  on  multiple  time  scales,  but  the  governing  effect  is  that  of  the  longest 
molecular  relaxation  time,  Tmax.  The  fluid  time  scale  is  estimated  as  the  reciprocal  of  the 
strain  rate  of  the  elongational  flow  [dU/dx]'1.  The  ratio  of  these  time  scales  is  called  the 
Weissenberg  number,  Wi,  and  is  used  as  a  measure  of  the  state  of  the  polymer  with 
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respect  to  the  coil-stretch  transition.  Often  the  elongational  time  scale  is  not  easily 
determined,  and  the  shear  time  scale  [dU/dy]'1  is  used  as  a  proxy. 

For  much  of  the  history  of  polymer  research,  only  bulk  properties  could  be  measured,  and 
the  behavior  of  polymers  as  a  function  of  Weissenberg  number  was  not  accurately 
known.  Wi  ~  1  was  taken  as  an  approximate  indicator  of  the  coil-stretch  transition. 
Beginning  in  1994,  techniques  for  single-molecule  observations  and  quantitative 
measurements  were  developed.  These  led  to  a  much  better  understanding  of  the 
dynamics  of  single  molecules,  and  to  an  appreciation  of  the  effects  of  different  types  of 
flows  on  the  coil-stretch  transition.  Key  papers  in  this  body  of  research  include  Perkins 
et  al,  1994a,  Perkins  et  al,  1994b,  Perkins  et  al,  1995,  Larson  et  al,  1997,  Perkins  et  al, 
1997,  Smith  &  Chu,  1998,  Smith  et  al,  1999,  Dua  &  Cherayil,  2000a,  Dua  &  Cherayil, 
2000b,  Hur  et  al,  2002,  Dua  &  Cherayil,  2003,  and  Schroeder  et  al,  2003. 

Most  pertinent  to  the  present  discussion  are  the  papers  of  Perkins  et  al,  1997  and  Smith  et 
al,  1999.  Comparison  between  those  two  papers  shows  that  there  is  a  very  significant 
difference  between  elongational  flow  and  shear  flow.  In  elongational  flow,  the  mean 
fractional  extension  of  the  polymer  molecule  increases  very  rapidly  with  Weissenberg 
number,  reaching  ~80%  mean  fractional  extension  for  Wi  ~  5,  then  slowly  approaching 
-95%  for  Wi  -  50.  In  shear  flow,  the  mean  fractional  extension  grows  much  more 
slowly  with  Wi  and  reaches  an  asymptotic  value  of  only  -  50%  for  large  Wi.  Figure  24 
below,  taken  from  Smith  et  al,  illustrates  the  aforementioned  point. 
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Figure  24:  |from  Smith,  et  al .,  Science  283,  1724-1727,  (1999)] 
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The  relevance  of  these  observations  for  slot-ejected  polymer  friction  reduction  can  be 
profound.  A  well-designed  ejection  slot  should  produce  primarily  elongational  flow,  in 
which  the  polymer  is  efficiently  extended  and  preconditioned  for  friction  reduction 
before  it  enters  the  external  flow.  Such  preconditioned  molecules  can  be  effective 
immediately.  On  the  other  hand,  if  the  ejection  slot  does  not  produce  elongational  flow,  a 
significant  fraction  of  the  polymer  will  enter  the  turbulent  boundary  layer  of  the  external 
flow  in  the  coiled  configuration.  As  Smith  et  al  showed,  the  shear  flow  of  the  boundary 
layer  is  much  less  efficient  at  uncoiling  the  polymer  molecules,  and  significant  delays 
may  ensue  before  friction  reduction  reaches  its  full  potential  (as  measured  by  Virk’s 
asymptote.) 

A  second  point  may  also  be  relevant.  Virk’s  asymptote  was  originally  derived  from 
empirical  observations  of  turbulent  shear  flow  in  pipes,  and  has  since  been  generalized  to 
turbulent  shear  flows  in  other  configurations.  It  represents  the  optimum  friction 
reduction  performance  that  can  be  achieved  in  shear  flow  with  a  mean  fractional 
extension  of  -50%.  When  the  mean  fractional  extension  of  the  polymer  molecules 
approaches  90%,  as  in  extensional  flow  at  large  Weissenberg  numbers,  the  friction 
reduction  performance  may  exceed  that  predicted  by  Virk’s  asymptote,  and  may  be 
described  by  a  separate  asymptote.  The  nature  of  that  separate  asymptote  is  yet  to  be 
defined,  but  one  can  conjecture  that  it  depends  on  the  Weissenberg  number  and  possibly 
on  the  Reynolds  number  as  well. 

A  preconditioned  polymer  solution  leaving  the  elongational  flow  environment  of  the 
ejection  slot  will  soon  interact  with  the  turbulent  external  flow.  As  a  result  of  mixing 
with  the  dominant  external  flow,  the  mean  fractional  extension  of  the  polymers  will  relax 
toward  50%,  consistent  with  shear  flow  conditions.  The  relaxation  time  for  this  process 
is  likely  to  be  significantly  longer  than  the  relaxation  time  for  a  single  molecule  (typically 
~  3  milliseconds),  but  the  exact  value  is  not  known.  The  relaxation  time  may  also  depend 
on  the  ejection  angle  of  the  slot,  which  would  influence  the  amount  of  time  spent  in  the 
viscous  sublayer  before  entering  the  turbulent  shear  flow  outside  the  viscous  sublayer. 

The  Weissenberg  number  for  the  ejection  flow  is  directly  proportional  to  q. 

Consequently,  the  mean  fractional  extension  of  the  ensemble  of  polymer  molecules  is  an 
increasing  non-linear  function  of  q,  consistent  with  Figure  24.  To  illustrate  this  point,  in 
a  simple  model  of  a  linearly  converging  channel  with  half-angle  a  and  exit  width  ws,  the 
maximum  rate  of  strain  at  the  exit  is 

dU/dxmax  =  2  tan(a)  uej/ws,  (4 1 ) 


and,  since 


q  =  n  Qs  =  Uej  ws, 


(42) 


the  Weissenberg  number  is 
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Wi  =  Tmax  (2  tan(a)/ws2)q. 


(43) 


Other  channel  geometries  would  preserve  the  linear  relationship  between  q  and  Wi,  albeit 
with  different  proportionality  factors. 

These  considerations  suggest  the  outlines  of  an  approximate  model  postulated  to  account 
for  the  near-slot  behavior  observed  in  the  Cortana  data. 

•  For  elongational  flow  there  exists  a  friction  reduction  asymptote  that 
exceeds  the  Virk  asymptote  for  shear  flow.  The  maximum  friction 
reduction  values  for  elongational  flow  are  presently  unknown,  but  are 
expected  to  be  a  non-linear  function  of  the  Weissenberg  number,  hence  of 
the  ejection  flow  rate. 

•  For  large  values  of  the  ejection  flow  rate,  local  friction  reduction  near  the 
slot  is  expected  to  attain  this  new  asymptote  almost  immediately. 

•  After  a  finite  period  of  mixing  with  the  external  flow,  the  behavior  of  the 
polymer  relaxes  from  the  new  asymptote  toward  the  Virk  asymptote  for 
shear  flow.  The  relaxation  time  is  equivalent  to  a  relaxation  distance, 
assuming  a  convection  velocity  that  is  some  fraction  of  the  free  stream 
velocity.  The  relaxation  time  is  presently  unknown,  but  can  be  estimated 
by  fitting  an  appropriate  model  to  the  near-slot  data. 

These  ideas  can  be  expressed  quantitatively.  The  relaxation  process  can  be  described  as: 

MAXDR  =  MAXDRvirk  ( 1  -  exp(-t/x))  +  MAXDREi0ng  cxp(-tA),  (44) 

where  t  is  the  relaxation  time  to  be  determined.  At  t  =  0,  MAXDR  =  MAXDREiong,  and 
as  t  — >oo,  MAXDR  —*  MAXDRvirk-  The  relationship  between  the  two  asymptotes  can  be 
derived  using  the  assumption  that  the  friction  coefficient  for  elongational  flow  is  less  than 
or  equal  to  the  Virk  asymptotic  friction  coefficient  by  a  factor  that  depends  on  the 
Weissenberg  number. 

Cfmong  =  F(Wi)  Cfyirk  0  <  F(Wi)  <  1  for  all  Wi.  (45) 

The  function  F(Wi)  =  F(q)  can  be  selected  to  fit  the  data.  These  assumptions  can  be 
combined  to  give 

MAXDR  =  MAXDRVirk  +  ( 1  -  F(q))(  1  -  MAXDRvirk)  exp(-t/x).  (46) 

Finally,  for  convection  downstream  from  the  slot, 

t  =  xs/U  =  qc0/pU2K,  (47) 
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MAXDR  =  MAXDRvirk  +  ( 1  -  F(q))(  1  -  MAXDRVirk)  exp(-qc0/pU2Kx).  (48) 


This  version  of  the  maximum  friction  reduction  asymptote  has  been  incorporated  in  the 
model,  with  the  value  of  x  and  the  specific  fonn  of  F(q)  to  be  determined  by  fitting  the 
Cortana  data  in  the  near-slot  region. 

Figure  25  shows  an  implementation  of  this  modification  to  the  model,  compared  with  the 
Cortana  data  taken  on  the  GTV  at  ARL  Penn  State.  The  experimental  data  are  the  same 
as  shown  in  the  earlier  Figure  24,  but  in  this  case  the  color  coding  is  different,  reflecting 
polymer  injection  rate  rather  than  free  stream  velocity.  Blue  is  now  for  q/Qs  =  5;  red  is 
for  q/Qs  =  10,  and  green  is  for  q/Qs  =  20.  The  model  has  been  adjusted  to  fit  the  data  by 
adjusting  x  and  defining  F(q).  The  fits  shown  in  the  figure  were  achieved  with  the 
relaxation  time  x  =  0.2  seconds  and  the  function  F(q)  defined  by  two  parameters: 

F(q)  =  exp(-((q/Qs)/15)4).  (49) 

This  functional  fonn  has  no  theoretical  significance,  but  gives  a  reasonable  fit  to  the  data. 
Note  that  the  relaxation  time  giving  a  reasonable  fit  is  1  to  2  orders  of  magnitude  greater 
than  the  molecular  relaxation  time,  consistent  with  observations  of  Pogrebnyak  & 
Ivanyuta  (1998)  for  high  concentrations  of  polymer. 

The  implication  of  the  data,  as  captured  in  the  model,  is  that  q/Qs  =  5  appears  to  provide 
very  little  elongation  of  the  WSR-301  polymer  in  solution,  whereas  q/Qs  =  20  appears  to 
provide  a  very  significant  amount  of  elongation. 
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Drag  Reduction 


K  Factor 

Figure  25:  WSR  301  Cortana  Fluidics  Ejector  on  GTV  at  ARL  Qs  Color  Coding 
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8  Model  Limitations  and  Additional  Research  Needs 


This  paper  has  shown  that  a  relatively  simple  phenomenological  model  can  predict  quite 
accurately  the  friction  reduction  performance  of  slot-ejected  polymer  solutions.  The 
model  contains  both  empirical  parameters  and  operational  parameters.  The  model  is 
particularly  sensitive  to  three  of  the  empirical  parameters;  the  power  law  exponent  [![! 
that  governs  the  approach  to  saturation  as  a  function  of  concentration  (for  low 
concentration),  the  slot  ejection  angle  expressed  through  the  parameter  a,  and  crou,  the 
value  of  concentration  at  which  friction  reduction  is  reduced  to  half  of  its  maximum  (for 
high  concentration).  Because  these  parameters  are  both  sensitive  and  poorly  understood, 
additional  research  should  focus  on  these  three  parameters. 

8.1  Improved  Estimates  of  pp 

The  nominal  estimate  pp  =  0.54  was  based  on  visual  data  extraction  from  a  figure  in  a 
published  paper.  The  data  set  was  very  limited  and  the  visual  extraction  was  probably 
not  very  accurate.  The  possible  range  for  the  parameter  was  estimated  as  0.4  <  PP  <  0.6. 

Estimates  of  pp  can  be  improved  by  examining  other  data  sets  for  which  both  the 
concentration  of  polymer  and  the  velocity  (or  Reynolds  number)  are  varied  over  a 
significant  range  and  are  accurately  measured.  Both  polymer  ocean  experiments  on 
smooth  plates  and  internal  flow  experiments  in  smooth  pipes  should  be  studied. 

The  underlying  assumptions  are  that  (1)  the  approach  to  saturation  of  friction  reduction 
follows  a  power  law  Co^/U*5,  where  pp  ^  8;  (2)  this  power  law  behavior  is  the  same  for 
polymer  ocean  experiments  and  pipe  experiments;  (3)  the  behavior  inferred  from  polymer 
ocean  experiments  and  pipe  experiments  can  be  generalized  to  slot-ejected  polymer 
experiments.  These  assumptions  have  some  support.  Assumption  (1)  appears  to  be  one 
of  the  necessary  conditions  for  the  separation  of  DR  vs.  K  curves  for  different  values  of 
Uoo  as  observed  by  Winkel  et  al. 

8.2  A  Model  for  a  vs.  <j) 

Many  past  experiments  have  used  ejection  angles  that  were  either  less  than  the  7°  or 
greater  than  the  20°  for  which  Vdovin  &  SmoTyakov  estimated  a.  Future  experiments 
will  also  exceed  those  limits.  A  reliable  means  for  extrapolating  a  would  be  very  helpful 
for  interpreting  such  experiments. 

More  important,  if  the  model  is  to  be  used  for  engineering  design,  it  must  be  capable  of 
accurate  prediction  of  friction  reduction  performance  for  different  (especially  very  small) 
ejection  angles.  The  goal  of  this  research  task  would  be  to  develop,  from  first  principles, 
a  fonnula  to  predict  a  as  a  function  of  cj).  It  is  likely  that  such  a  fonnula  would  require  a 
more  comprehensive  approach  to  the  diffusion  of  polymer  in  the  near-wall  region,  and 
might  lead  to  deeper  insights  about  the  diffusion  length,  L,  introduced  by  Vdovin  & 
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Smol’yakov.  In  particular,  it  might  clarify  why  L  behaves  differently  for  small  and  large 
ejection  angles  and  why  L  appears  to  approach  an  asymptote  Lmax under  some  conditions. 

As  noted  above,  a  critical  and  fundamental  issue  is  the  observation  of  Vdovin  & 
Smol’yakov  that,  for  small  ejection  angles,  the  behavior  of  the  polymer  is  governed  by 
inner  variables  alone,  not  by  outer  variables  such  as  the  free  stream  velocity. 
Understanding  in  more  detail  the  basic  physics  behind  this  observation  will  be  a 
significant  advance. 

8.3  Reliable  Estimate  for  croii 

The  value  of  concentration  for  which  friction  reduction  performance  drops  by  M>  needs  to 
be  estimated  more  accurately,  and  the  physics  that  produces  this  high-concentration  roll¬ 
off  needs  to  be  better  understood.  The  empirical  data  of  Wu,  Fruman  &  Tulin  and  of 
Povkh  et  al  are  qualitatively  similar,  but  do  not  completely  agree  quantitatively.  The 
functional  form  of  the  roll-off  used  in  the  model,  croii/(c  +  cron),  is  only  a  rough 
approximation  that  captures  the  general  behavior.  In  addition,  Povkh  et  al  show  that  the 
roll-off  behavior  depends  on  other  parameters  such  as  molecular  weight  that  are  not 
included  in  the  present  model. 
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9  Summary 


The  goal  of  this  work  was  to  develop  a  model  that  accurately  captures  the  drag-reduction 
behavior  of  slot-ejected  polymer  in  external  flow  and  that  can  be  used  as  a  measure  of 
effectiveness  (MOE)  for  engineering  purposes. 

An  empirical  model  has  been  developed  by  combining  partial  models  from  the  work  of 
earlier  researchers.  Models  of  polymer  wall  concentration  as  a  function  of  downstream 
distance  from  the  slot  were  taken  from  Vdovin  &  SmoTyakov  (1978,  1982).  A  model  of 
the  dependence  of  friction  reduction  on  wall  concentration  was  inferred  from  data 
reported  by  Petrie  et  al  (2003).  Models  for  maximum  achievable  friction  reduction  as  a 
function  of  Reynolds  number  were  taken  from  Virk  et  al  (1970)  and  from  White  (1991). 
Data  on  the  decrease  of  friction  reduction  perfonnance  at  high  polymer  concentration 
were  taken  from  Wu,  Fruman  &  Tulin  (1977)  and  from  Povkh  et  al  (1979).  These  model 
components  were  combined,  using  simple  subsidiary  calculations,  to  create  a  complete 
model  of  friction  reduction  as  a  function  of  the  dimensionless  parameter  K. 

The  complete  model  contains  several  empirical  parameters  from  earlier  published  works 
and  several  operating  parameters  that  define  an  experimental  configuration.  Each  of  the 
parameters  was  tested  for  its  impact  on  the  model.  It  was  found  that  the  model  was  very 
sensitive  to  three  parameters  in  particular:  the  power  law  exponent  that  governs  the 
approach  to  saturation  as  a  function  of  polymer  concentration,  the  ejection  angle  of  the 
slot  that  extrudes  the  polymer  into  the  flow,  and  the  roll-off  concentration  at  which 
perfonnance  decreases  for  high  polymer  concentrations.  The  power  law  parameter 
controls  the  slope  of  the  DR  vs.  K  curve  in  the  end  diffusion  region.  The  ejection  angle 
controls  the  extent  of  the  intermediate  diffusion  region,  and,  consequently,  has  a  major 
impact  on  overall  performance.  The  high-concentration  roll-off  parameter  controls 
behavior  close  to  the  ejection  slot  where  initial  polymer  concentration  is  high. 

For  large  values  of  K  (regions  near  the  slot)  performance  is  generally  controlled  by  the 
Virk  asymptote  and  the  high  concentration  roll-off.  For  small  values  of  K,  perfonnance 
is  controlled  by  the  end  diffusion  zone.  Between  these  two  limits,  performance  may  be 
controlled  by  the  intennediate  diffusion  zone.  However,  for  large  ejection  angles,  the 
width  of  the  intennediate  zone  becomes  insignificant,  and  performance  transitions 
directly  from  the  Virk  asymptote  region  to  the  diffusion  far  zone  region.  The 
intennediate  diffusion  zone  is  important  only  for  small  ejection  angles. 

The  model  was  fine-tuned  and  compared  with  data  from  published  experiments. 
Excellent  agreement  with  data  from  Winkel  et  al  was  demonstrated.  In  particular, 
differential  perfonnance  due  to  different  free  stream  velocities  and  to  different  polymer 
types  was  predicted  accurately.  Models  of  friction  reduction  performance  at  shallow 
ejection  angles  produced  good  agreement  with  data  taken  at  ARE,  Penn  State  for  Cortana 
Corporation.  An  adjunct  model  was  developed  to  explain  anomalous  behavior  in  the 
near-slot  region. 
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Limitations  of  the  model  were  discussed.  Three  areas  of  high  priority  research  were 
identified:  (1)  more  accurate  estimation  of  the  power  law  parameter  for  the  low- 
concentration  approach  to  saturation  and  the  roll-off  parameter  for  high  concentration;  (2) 
a  first  principles  model  of  the  effect  of  changing  the  ejection  angle;  (3)  additional 
modeling  of  near-slot  behavior  of  dilute  polymer  solutions,  considering  the  effects  of 
concentration  and  molecular  weight. 

The  model  provides  a  useful  MOE  tool  for  evaluation  of  field  experiments  and  for  design 
of  slot  ejector  systems.  The  local  DR  vs.  K  curves  can  be  inverted  to  predict  local  DR  vs. 
downstream  distance.  That  function  can  then  be  integrated  to  predict  total  friction 
reduction  as  a  function  of  polymer  consumption  for  various  slot  designs.  Performance 
for  single  slots  or  multiple  slots  can  be  evaluated  and  compared.  In  this  way,  an  optimum 
polymer  ejection  system  can  be  developed,  its  performance  can  be  evaluated,  and  its 
cost/effectiveness  vis-a-vis  reduced  fuel  consumption  can  be  determined. 
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11  Glossary 


Alpha-  Alpha  (a)  is  a  parameter  that  depends  on  angle  at  which  the  ejectant  is  introduced  into 
the  free  stream  flow. 


Beta  -  Beta  (P)  is  a  parameter  that  depends  on  the  ejection  angle  of  the  slot,  p  =  1  -  a  to  ensure 
that  cw/co  =  1/e  when  x  =  L,  defined  as  the  diffusion  length. 

Boundary  Layer  -  The  boundary  layer  (BL)  is  defined  as  the  region  in  an  external  flow  where  the 
velocity  varies  from  zero  at  the  wall  to  approximately  99  percent  of  the  free  stream 
velocity.  It  is  the  region  where  the  effects  of  fluid  viscosity  (r|)  are  present.1 2 

Buffer  Layer  -  The  region  where  both  viscous  and  Reynold’s  stresses  are  comparable  (also 
referred  to  as  transition  layer). 


Correction  Factor  -  The  correction  factor  (F)  is  a  factor  to  correct  for  the  maximum  diffusion 

kqc:0/ 

length  (Zmax)  and  is  given  by  the  equation:  F  =  — / 

K<ico 

^  _  ePU«Lm ax 


Delta  -  In  this  context,  the  parameter  delta  (8)  is  the  power  law  exponent  that  describes  the 
contribution  of  the  velocity  in  the  approach  to  saturation  as  a  function  of  wall 
concentration. 

Density  -  The  mass  per  unit  volume  is  represented  by  the  Greek  letter  rho  (p).  For  seawater,  its 
mass  density  at  28.3  °C,  p  =  1022.2  kg/m3. 

Diffusion  Length  -  The  diffusion  length  ( L )  as  defined  by  Vdovin  and  Smol’yakov  is  the 

downstream  distance  at  which  the  wall  concentration  has  dropped  to  1/e  of  its  initial 
value  (where  e  is  the  base  of  the  natural  logarithm). 

Double  Beta  -  The  parameter  double  beta  (pp)  is  the  power  law  exponent  of  the  concentration  in 
the  approach  to  saturation  from  lower  concentrations.  Numerically,  it  is  a  constant 
having  a  value  of  less  than  one. 

Downstream  Distance  -  The  downstream  distance,  denoted  by  the  tenn  (x)  is  relative  to  the 
virtual  origin  such  that  x  =  x0  +  xs  when  x0  is  the  coordinate  of  the  slot,  and  xs  is  the 
downstream  distance  from  the  slot. 


1  Gillmer,  Thomas  C.  and  Bruce  Johnson,  Introduction  to  Naval  Architecture,  Naval  Institute  Press, 
Annapolis,  MD,  1982,  pg  216. 

2  Tennekes,  H.  and  J.L.  Lumley,  A  First  Course  in  Turbulence,  the  MIT  Press,  Cambridge,  MA  and 

London,  England,  16th  printing  1997,  pg  161. 
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Friction  reduction  -  For  this  paper,  friction  reduction  will  refer  to  a  reduction  of  the  skin  friction 
drag  which  is  caused  by  shear  stresses  at  the  surface.  Again,  for  this  paper,  the  reduction 
in  skin  friction  drag  is  attained  via  polymer  ejection. 

Dynamic  Viscosity  -  “A  measure  of  the  resistance  to  flow  of  a  fluid  under  an  applied 
force”  denoted  by  the  Greek  symbol  mu  (|ii  =  vp). 

Ejection  Angle  -  The  angle,  usually  measured  in  degrees,  at  which  the  ejected  solute  is 

introduced  into  the  free  stream  via  the  ejector.  This  angle  (<|>)  is  typically  measured  as  the 
angle  between  the  upstream  free  stream  surface  and  the  upstream  ejector  surface. 

Ejection  Velocity  -  The  velocity  at  which  the  solute  is  ejected  from  the  ejection  orifice  into  the 
free  stream  (uej). 

Far  Zone  -  As  the  turbulent  boundary  layer  (TBL)  thickens  downstream  from  the  point  of 

ejection,  diffusion  continues  to  reduce  the  concentration  of  solute  near  the  wall,  but  the 
rate  is  algebraic  rather  than  exponential.  On  the  K-Plot,  the  Far  Zone  is  also  the  region  of 
changing  slope  to  the  left  of  the  constant-slope  region,  approaching  zero  drag  as  K 
decreases.  The  Far  Zone  is  also  sometimes  referred  to  as  the  Final  Zone  or  End  Zone. 

Gamma  -  The  term  gamma  (y)  is  a  constant  that  depends  on  the  flow  velocity.  The  parameter  y 
appears  in  the  exponent  that  governs  the  approach  to  saturation  as  a  function  of  wall 
concentration.  In  general,  increasing  y  increases  the  apparent  perfonnance  of  the 
polymer. 

Initial  Concentration  -  Typically  refers  to  the  concentration  of  the  polymer  powder  in  the  fluid 

being  ejected  on  a  weight  basis.  The  initial  concentration  is  often  denoted  by  the  symbols 
Co  or  Ci  and  is  usually  described  in  weight  parts  per  million  (wppm). 

Initial  Zone  -  Located  very  close  to  the  slot,  where  the  polymer  is  confined  to  the  near-wall 

region  and  the  wall  concentration  remains  close  to  the  initial  concentration.  Typically, 
the  Initial  Zone  is  very  short. 

Intermediate  Zone  -  In  this  zone,  the  diffusion  process  transports  the  solute  (in  this  case,  the 

polymer  additive)  away  from  the  wall  until  it  occupies  much  of  the  turbulent  boundary 
layer.  The  concentration  of  the  solute  at  the  wall  decreases  exponentially  in  this  region  or 
zone.  The  Intermediate  Zone  is  also  sometimes  referred  to  as  the  Transition  Zone. 

ko  -  Performance  factor  for  polymer  as  derived  from  Vdovin  and  Smol’yakov’s  seminal  work. 
They  established  a  method  for  comparing  DR  results  that  used  different  types  of 
polymers.  The  reference  value  for  which  ko  is  equivalent  to  one  (1)  is  Dow  Chemical’s 
PEO  (PolyEthylene  Oxide)  WSR-301. 


3  http:// dictionary. reference. com/browse/dvnamic%2QViscosity 


52 


k  -  This  term  is  derived  from  the  reference  ko.  Polymers  that  have  twice  as  much  ability  to 

provide  DR  as  WSR-301  would  have  a  k-value  of  two  (2)  and  polymers  that  were  only  Vi 
as  efficient  as  WSR-301  would  have  a  k-value  of  0.5. 


K-Curve  -  A  curve  plotted  on  a  K-factor  graph  which  usually  has  the  K-factor  plotted 

logarithmically  on  the  x-axis  (abscissa)  and  %  DR  (specifically,  this  is  percentage  friction 
reduction,  not  percentage  drag  reduction)  linearly  plotted  on  the  y-axis  (ordinate). 


K-Factor  -  A  value  based  on  the  work  by  Vdovin  &  SmoFyakov 


K  -  Factor  = 


Q*Cj 

P*UxxXle 


2  •  *3 

where  Q  =  the  ejection  mass  flow  rate  (crrr/scc),  Q  =  the  ejectant  concentration  (g/cm  ), 
p  =  the  density  of  the  fluid  medium  (g/cm3),  U  «,  =  the  velocity  of  the  external  flow 
(cm/sec),  and  Xte  =  the  distance  from  the  ejection  slot  to  the  trailing  edge  of  the  drag 
balance  (cm). 


Kinematic  Viscosity  -  This  is  a  coefficient  which  describes  the  diffusion  of  momentum.4 5  The 
kinematic  viscosity  (v)  is  calculated  taking  the  dynamic  viscosity  (viscosity  coefficient 
that  determines  the  dynamics  of  incompressible  Newtonian  fluid)  and  dividing  by  the 


density  of  Newtonian  fluid 


5 


Laminar  Flow  -  “Sometimes  known  as  streamline  flow,  occurs  when  a  fluid  flows  in  parallel 

layers,  with  no  disruption  between  the  layers.  In  fluid  dynamics,  laminar  flow  is  a  flow 
regime  characterized  by  high  momentum  diffusion,  low  momentum  convection,  pressure 
and  velocity  independent  from  time.  It  is  the  opposite  of  turbulent  flow.  In  nonscientific 
tenns  laminar  flow  is  smooth,  while  turbulent  flow  is  rough.”6  “In  laminar  flow,  the  fluid 
moves  in  smooth  layers  or  lamina.  There  is  relatively  little  mixing  and  consequently  the 
velocity  gradients  are  small  and  shear  stresses  are  low.  The  thickness  of  the  laminar 
boundary  layer  increases  with  distance  from  the  start  of  the  boundary  layer  and  decreases 
with  Reynolds  number.”7 

Length  -  The  measure  of  the  greatest  dimension  of  a  plane  or  solid  figure  often  denoted  by  the 
symbol  (L)  and  is  usually  in  units  of  ft,  meters,  etc. 

Length  of  the  Plate  -  The  length  of  the  plate  (Lp)  refers  to  a  surface  with  an  infinite  extent  in  the 
direction  perpendicular  to  the  flow.  The  surface  may  be  either  a  flat  plate  or  a  body  of 
revolution  having  length,  Lp. 

Maximum  Friction  Reduction  -  For  a  specific  case  (i.e.  ejector  type,  solute  flow  rate,  external 
velocity, ,  this  is  the  maximum  amount  of  friction  reduction  that  was  attained. 


4  http://scienceworld.wolfram.com/phvsics/KinematicViscosity.html 

5  http://en.wikipedia.org/wiki/Kinematic  Viscosity 

6  http://en.wikipedia.org/wiki/Laminar  flow 

7  http://www.desktopaero.com/appliedaero/blavers/blayers.html 
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Reynolds  Number  -  The  Reynolds  number  (Re),  refers  to  “the  ratio  of  inertial  forces  ( vsp )  to 
viscous  forces  (ju/L)  and  consequently  it  quantifies  the  relative  importance  of  these  two 

o 

types  of  forces  for  given  flow  conditions.” 

Roll-off  Parameter  -  An  empirical  parameter  for  high-concentration  effects  that  denotes  the 

movement  of  solute  from  an  area  of  higher  concentration  near  the  wall  to  an  area  of  lower 
concentration  further  out  into  the  boundary  layer  (croii). 

Saturation  -  The  non-linear  approach  of  a  dependent  variable  to  a  fixed  asymptote  as  an 

independent  variable  increases  linearly.  In  this  paper,  the  terminology  as  applied  to  the 
approach  of  friction  reduction  to  its  maximum  value  as  polymer  concentration  increases. 

Slot  Width  -  The  slot  width  refers  to  the  ejector  orifice  as  measured  from  the  upstream  side  of  the 
ejector  opening  to  the  downstream  side  of  the  ejector  (ws).  Ejectors  widths  are  typically 
on  the  order  of  1  mm  across. 

Tau  -  The  term  tau  (t)  refers  to  the  relaxation  time  of  the  polymer  molecule. 

Toms  Effect  -  In  1949,  the  British  researcher  B.A.  Toms  discovered  that  dissolving  small 

quantities  of  heavy,  long-chain  molecular  polymers  in  a  solution  reduced  the  drag  in  a 
turbulent  flow  by  upwards  of  70  percent. 

Turbulent  Flow  -  “In  fluid  dynamics,  turbulence  or  turbulent  flow  is  a  flow  regime  characterized 
by  chaotic,  stochastic  property  changes.  This  includes  low  momentum  diffusion,  high 
momentum  convection,  and  rapid  variation  of  pressure  and  velocity  in  space  and  time.”9 

Turbulent  Boundary  Layer  -  “The  turbulent  boundary  layer  (TBL)  flow  is  characterized  by 
unsteady  mixing  due  to  eddies  at  many  scales.  The  result  is  higher  shear  stress  at  the 
wall,  a  "fuller"  velocity  profile,  and  a  greater  boundary  layer  thickness.  The  wall  shear 
stress  is  higher  because  the  velocity  gradient  near  the  wall  is  greater.  This  is  because  of 
the  more  effective  mixing  associated  with  turbulent  flow.  However,  the  lower  velocity 
fluid  is  also  transported  outward  with  the  result  that  the  distance  to  the  edge  of  the  layer 
is  larger.”10 

Velocity  -  A  vector  quantity  whose  magnitude  is  a  body's  speed  and  whose  direction  is  the 

body's  direction  of  motion.* 11  Referred  to  as  (Uref)  with  the  units  kts,  ft/sec,  or  m/sec  and 
can  also  refer  to  the  free  stream  velocity  (Uqo). 

Viscous  Sublayer  -  This  term  refers  to  a  zone  of  unsteady  flow  where  velocity  fluctuations  do  not 

12 

contribute  much  to  the  total  stress  due  to  the  overwhelming  effects  of  viscosity. 


8  http://en.wikipedia.org/wiki/Reynolds  number 

9  http://en.wikipedia.org/wiki/Turbulent  flow 

10  http://www.desktopaero.com/appliedaero/blavers/blayers.html 

11  http://dictionarv.reference.com/browse/velocity 

12  Tennekes,  H.  and  J.L.  Lumley,  A  First  Course  in  Turbulence,  the  MIT  Press,  Cambridge,  MA  and 
London,  England,  16th  printing  1997,  pg  160. 


54 


Virtual  Origin  -  Denoted  as  the  position  (x)  along  the  plate  at  which  the  turbulent  boundary  layer 
begins.  The  virtual  origin  is  assigned  the  coordinate  x  =  0,  and  downstream  positions  are 
referenced  to  that  origin. 

Volume  Ejection  Rate  per  Unit  Span  -  This  term  is  the  volume  discharge  per  unit  span  width 

flow  rate  of  the  ejected  additive  (solute)  into  the  surrounding  fluid  medium  (typically  sea 
water  in  an  at-sea  experiment)  and  is  denoted  in  units  of  quanta  of  viscous  sublayer 

volume  flow  rates  [viscous  sublayer  flow  rate  0<y+U  1.6  and  is  =  67. 3v  in  ft  /sec  or 

TIT. 

m  /sec].  This  term  for  the  volume  ejection  rate  per  unit  span  is  often  denoted  by  either 
the  symbol  Qs  or  q. 

Weissenberg  Number  -  “A  dimensionless  number  used  in  the  study  of  viscoelastic  flows.  It  is 
named  after  Karl  Weissenberg.  The  dimensionless  number  is  the  ratio  of  the  relaxation 
time  of  the  fluid  and  a  specific  process  time.  For  instance,  in  simple  steady  shear,  the 
Weissenberg  number,  often  abbreviated  as  Wi  or  We,  is  defined  as  the  shear  rate  times 
the  relaxation  time  (Wi  =  yX).  Since  this  number  is  obtained  from  scaling  the  evolution 
of  the  stress,  it  contains  choices  for  the  shear  or  elongation  rate,  and  the  length-scale. 
Therefore  the  exact  definition  of  all  non  dimensional  numbers  should  be  given  as  well  as 
the  number  itself.”14 


13  Fontaine,  A.A.,  FI.L.  Petrie,  and  T.A.  Brungart,  “Velocity  Profile  Statistics  in  a  Turbulent  Boundary 
Layer  with  Slot-Injected  Polymer”  Journal  of  Fluid  Mechanics  (1992),  vol.  238,  pp  435-466. 

14  http://en.wikipedia.org/wiki/Weissenberg  number 
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